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A dense vortex lattice in a rotating dilute Bose-Einstein condensate is studied with the Thomas- 
Fermi approximation. The upper critical angular velocity fi C 2 occurs when the intervortex separation 
b becomes comparable with the vortex core radius £. For a radial harmonic trap, the loss of con- 
finement as — » lu± implies a singular behavior. In contrast, an additional radial quartic potential 
provides a simple model for which Q C 2 is readily determined. Unlike the case of a type-II super- 
conductor at fixed temperature, the onset of fl C 2 arises not only from decreasing b but also from 
increasing £ caused by the vanishing of the chemical potential as Q — > Q C 2- 
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I. INTRODUCTION 

The recent observations [^-^ of vortex arrays and vortex lattices in rotating dilute trapped Bose-Einstein conden- 
sates (BEC) raise the question of analogies with type-II superconductors. What limits the number of vortices? Is 
there an upper critical rotation speed fl C 2 similar to the upper critical field H C 2 [^) familiar from type-II supercon- 
ductors? In the superconducting case, the field H C 2 occurs when the distance b ~ ^JQq/ttB between vortices becomes 
comparable with the vortex core size £ ~ 10-100 run in typical conventional superconducting alloys. The measured 
flux quantum $0 = h/2e then implies H C 2 ~ ^o/ 71 "^ 2 ~ 0.1-10 T (in this limit, the distinction between H and B 
becomes negligible). For rotating superfluid 4 He where £ ~ a few A and b ~ yJhjM^l 0], essentially the same 
condition yields the unattainably large value Q C 2 ~ h/M^ 2 ~ 10 15 rad/s. 

A similar criterion also applies to rotating Bose condensates. In this case, the vortex core radius £ = h/y/2M/i ~ 
0.1 fjxa is macroscopic, where \x is the chemical potential, suggesting that experimental study of f2 c2 might well be 
possible. As seen below, the situation in rotating dilute condensates is even more favorable, because the chemical 
potential decreases at large rotation speeds, increasing the vortex-core size and thus reducing Q, C 2- 

For a trap with a radial harmonic potential V±(r±) = hMuJ\r\, the behavior becomes singular when ft — > u± 
because the outward centrifugal force counteracts the inward force from the harmonic trap 0-0] • Specifically, the 

effective radial trap frequency {u>\ — fl 2 ) produces an observable centrifugal distortion of the condensate Jl(]-[lJ] 
for currently attainable values of fi/cjj_ 1. As seen below, this behavior means that uj± effectively acts like f2 c2 
for pure harmonic radial confinement, and direct experimental study of the limiting behavior for f2 — * uj±_ would be 
difficult. Thus, it is convenient to consider an additional stiffer radial potential, which eliminates the singularity when 
il = uj_\_. In this case, £l c2 exceeds uj±, and the limit Q — > £l C 2 is relatively smooth. The analysis is especially simple 
for a quartic radial potential Vi(r±) = \kr 4 Ll and this example is analyzed in detail. 

Section II introduces the basic Gross-Pitaevskii free energy and condensate wave function ^ = |^|e for a rotating 
vortex lattice in a trapped BEC. In equilibrium, the vortex lattice must experience self-consistent solid-body rotation, 
which means that neither the phase S nor the superfluid velocity v s = (hfM)'VS can be spatially periodic in the 
laboratory frame. The transformation to the frame rotating with angular velocity f2 yields a solid-body velocity 
v s b =!)xr that cancels the overall rotation. Thus, the resulting relative velocity v s — v s b is indeed spatially periodic 
for an unbounded vortex lattice. This feature facilitates a simple description of the vortex lattice in a large trapped 
BEC. The case of a trap with harmonic radial confinement is studied in Sec. Ill, and the more interesting case of an 
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additional quartic radial confining potential is studied in Sec. IV 



II. GENERAL FORMALISM 



It is convenient to start from the Gross-Pitaevskii (GP) JT^,[l4| free energy in the rotating frame 



F = dV 



2M 
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where V±(r±) is the transverse radial confining potential and g = Airafi 2 /M relates the interparticle coupling constant 
to the s-wave scattering length a (here assumed positive). The representation \& = | \I/ 1 e 5 " 5, emphasizes the hydrody- 
namic aspects of the behavior, with condensate density n = |\1/| 2 and superfluid velocity v s = h'VS/M. The quantity 
— vff* ft ■ r x p*f? can be written as ^iMl d \^>\ 2 /d<p — M v s b • v s |^| 2 , and the first term makes no contribution to the 
spatial integral. Straightforward manipulations of Eq. (Q) yield 



F = dV 
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where V ccn t(r±) = —\M$}r\ is an effective repulsive centrifugal potential. 

The simplest variational model is to take the phase S as that for a classical unbounded vortex lattice (with vortices 
aligned along z), which has been analyzed in detail by Tkachenko |15| (as a model for a vortex lattice in rotating 
superfluid 4 He). In this case, the total superfluid velocity includes the divergent locally axisymmetric circulating flow 
near each vortex core. The resulting singular kinetic energy in Eq. (g) is cut off by a self-consistent core structure of 
characteristic dimension £ = h/y/2Mjl, where /x is the chemical potential fl3| , |l4| , |l6f . 

This classical hydrodynamic phase S cannot be a spatially periodic function in the xy plane because the vortex 
array induces a self-consistent rotation. Instead, this phase obeys definite quasiperiodic conditions |]l5| , [i"7f that depend 
on the details of the lattice structure. For the same reasons, v s cx VS* is also not spatially periodic. In contrast, 
the quantity v s — v s b (namely the superfluid velocity as seen in the rotating frame) is spatially periodic, because 
the subtracted term v s b cancels the effect of the quasiperiodic terms in VS, considerably simplifying the subsequent 
analysis. 

According to Feynman's picture of a rotating superfluid Jl8| , the vortices have a uniform areal density n v = 2MQ/h, 
ensuring that the mean vorticity of the vortex lattice mimics that of solid-body rotatio n V x v s b = 2fi. Thus the area 
per vortex n" 1 = irb 2 can be taken to define an intcrvortex separation b = uK/MO,] this characteristic length sets 
the scale of the vortex lattice. For the present purposes, the detailed lattice structure is unimportant (for example, 
the free energies of the triangular and square configurations have the same logarithmic contributions oc ln(6/£) and 
differ only in the additive constants p9[); it will be convenient to use a Wigner-Seitz approximation in which each 
polygonal unit cell is replaced by an equivalent circular cell of radius b. Evidently, the ratio £ 2 /& 2 characterizes the 
fractional volume occupied by the "normal" vortex cores. 

For any reasonable transverse confining potential V±(r±), the single-particle ground state will have some charac- 
teristic transverse dimension d± , alongwith the characteristic axial dimension d z — yjTi/Muj z that is set by the axial 
harmonic confining potential in Eqs. (|l|) and (^). Let Nq be the number of atoms in the condensate at low temper- 
ature. In the Thomas- Fermi (TF) limit ^0|, the condensate density is taken as locally constant, which holds when 
the dimensionless interaction parameter Noa/d± is large. In this case, the repulsive Hartree (mean-field) interactions 
expand the condensate relative to its noninteracting transverse and axial dimensions d± and d z . In particular, the 
radial and axial condensate radii R± and R z are simply the classical turning points for a particle of energy //. Thus, 
R z = y // 2fj,/Mu> z depends only on ^t, but R± also depends on fl because of the repulsive centrifugal potential V CC nt (rx ) 
in Eq. (§. 

For a rotating condensate with chemical potential [i 1 the GP equation derived from Eq. (^) implies the corresponding 
TF density 
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obtained by omitting the gradient of the condensate density. If the vortex cores do not overlap significantly (so that 
£ 2 /6 2 <C 1), the condensate density can be approximated by a product ^(r^z)! 2 w \^tf(t±, z)\ 2 u 2 (r±) of the TF 
density 



\^TF(rx,z)\ 2 = ^ 
9 

in the absence of the vortex lattice, and a factor 



Vj_(rs_) tt 2 r 2 ± Mlu 2 z z 2 
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that is spatially periodic in the transverse plane and must be cut off near the vortex core to ensure that u 2 remains 
positive. In any given unit cell, it has the local form 
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where a circular Wigner-Seitz cell has been used with £ < r j_ < b. In effect, the resulting density is that of a vortex-free 
condensate with narrow holes along the axes of the vortex lattice . 

The chemical potential is determined by the general condition that No — J dV \ ^(r±, z)\ 2 . Here the factorized form 
yields 



No- J dV\^TF(r±,z)\ 2 u 2 (r ± ) 
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For a dense vortex lattice, the intervortex separation b is small compared to the TF transverse radius R± , so that the 
first factor varies slowly on the length scale b. Thus, u 2 can be averaged over any single unit cell, and the Wigner-Seitz 
approximation JT^ ] in Eq. (^|) yields 



(u 2 ) 
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as anticipated, (u 2 ) ss 1 for £ <C b. In the limit of a large condensate with many vortices, Eq. (Q) can therefore be 
approximated by 
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III. RADIAL HARMONIC CONFINING POTENTIAL 



The simplest situation is the harmonic radial potential V±(r±) = in which case Eq. (@) implies the TF 

condensate radius 



Familiar manipulations yield 



^W mk%) - (10) 
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where A = w z /cJj_ is the axial anisotropy parameter. 

Except for a very narrow region £l/u>± < 1 close to the singular limit, the vortex cores occupy negligible volume, 
so that (u 2 ) ss 1, In this case, the condensate particle number No remains essentially constant, and Eq. ([n]) shows 
how the chemical potential decreases with increasing rotation speed Q = Q/uj±. In this same large interval, Eq. (|Io| ) 
implies that 

Rl(Q) _ /157V aA\ 2/5 A n 2\-3/5 
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Equivalently, 

i?±(0) 

as seen in recent experiments [ flO|| . 

The fraction of depleted condensate occupied by the vortex cores is given quite generally by £ 2 jb 2 = 7if2/2/x, and 
use of Eq. (O) gives the specific result 



^^-^-0 2 )- 3/1 °, (13) 
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where 



b2 ^ 2 (i-n 2 ) 275 



d± «1 (15) 



15iV aA R\{Q) 

is small in the TF limit. Equation ( |l4| ) shows that the ratio S, 2 /b 2 increases linearly with Q but remains small (of 
order e 2 / 5 ) until 1 — Q < e/8v% when the ratio grows rapidly toward 1. If f2 C 2 is defined like H C 2 for a type-II 
superconductor, f2 C 2 corresponds to the limit £ ~ 6. As a result, a radial harmonic trap potential leads to singular 
behavior when fl — > c<j_l |^] . Specifically, the disappearance of the condensate associated with the sudden expansion 
of the "normal" vortex cores presumably implies a phase transition, which occurs essentially simultaneously with the 
loss of confinement as Q — ► 1 |Q-^,^2) . 

This straightforward analysis indicates that any study of a dense vortex lattice near fi C 2 for a radial harmonic 
trap will encounter significant difficulties arising from the softening of the effective trap potential. As a result, it is 
interesting to consider a stiffer radial potential, and the next section treats one possible example. 



IV. RADIAL QUARTIC CONFINING POTENTIAL 

As a specific example, let V±(r±) include a quartic radial trap potential Vi(r±) = \kr\ in addition to the quadratic 
potential \Mlo 2 j 2 _\ Lundh |2^] has recently examined this and other power-law potentials in a theoretical study of 
the formation of a multiply quantized vortex in a rotating condensate. For a two-dimensional ideal gas confined 
with pure Vi(r±), the balance between the ground-state kinetic energy ~ Ti 2 /Md 2 and ground-state potential energy 
~ kd\ implies a ground-state mean radius du ~ (h 2 /Mk) 1 ^ 6 and ground-state energy E4 = huj^ ~ (Ti k/M 2 ) 1 / 3 . For 
definiteness, it is convenient to ignore numerical constants of order unity and take 



Ea = Tiloa 



Md 2 ~ kd ' 4 ' 



(16) 



as the relevant dimensional quantities. 

For a rotating condensate in this combined radial trap, Eq. ([)]) becomes 
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where u 



and R\ is the turning point where the integrand vanishes. The presence of the quartic potential means 
that the external rotation speed £1 can now exceed u>±. In this regime (11 > 1), the particle density actually attains 
a local minimum on the axis of symmetry, but this effect is probably difficult to detect. 
The integral in Eq. (ITTn can be expressed in terms of a dimensionless parameter 
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A detailed analysis yields the final form 



du 



m{uo\ - n 2 ) 

2[L 



u u 

4 M 



3/2 



2 VV<»>- 



(19) 



where 



3tt 
16 



3?7 3 



5// 
~8' 



(20) 



It has the limiting forms 
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(5t7) _1 for r) > 1, 
f(v)~{ TE«-V for H<1, (21) 



§7r?7 4 for 77 < -1. 



Equations ( |l7| ) and ([19]) can be combined with R z = ^2/i/Mu 2 to express the chemical potential in terms of the 
parameter 77 
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where d±/di is a dimensionless measure of the relative strength of the quartic potential (note that d±/di — > in 
the limit that the quartic coupling constant k becomes small). In addition, Eq. (|l^) can be rewritten to express the 
rotation speed as a function of /1 and r\ 

tf^-l-^JjLl^ (23) 



uj" 2 ^ 1 y fixjj 

Note that 77 is large and positive for a nonrotating condensate (Q = 0), r\ vanishes for fi = 1, and 77 becomes negative for 
£7 > 1 (this limit can occur only for nonzero positive d±/di). These two equations provide a parametric representation 
of the dependence of the dimensionless chemical potential JL = fi/huj± on the rotation speed, generalizing the result 
in Eq. (|llj) for a pure harmonic potential to include the effect of an additional quartic potential. 

It is not difficult to see from Eqs. ( pl| ) and ( p2] ) that /1 vanishes like 77" 2 as 77 — * —00, so that the right-hand side of 
Eq. (^) remains finite in the same limit. Thus the chemical potential vanishes when the rotation frequency attains 
its maximum value. This maximum rotation speed can be identified as the upper critical value f2 C 2 because the vortex 
core size £ = Ti/ ^/2M /1 diverges as /1 — + 0. A detailed analysis (approximating (u 2 ) « 1) shows that 

\ 7T di UJ4 I 

where the positive quantity fl 2 2 — cjj_ has been expressed solely in terms of the parameters associated with the 
quartic term in the confining potential. Figure [l] shows the dimensionless chemical potential /TZ as a function of the 
dimensionless rotation speed f2 for two illustrative cases (note that ~p remains positive for O = 1). 
In the presence of the additional quartic confining potential, the TF radius R± is given by 



R ± _ o ( <** x;: 



d\ V huJ ^- \ d ^J \/l + rf + 77 



(25) 



A combination with Eq. (E3h provides a parametric relation for R±(Q). Figure shows the dimensionless ratio R\jd 
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as a function of the dimensionless rotation speed O (note that Rj_/dj_ remains finite as Q — > fl C 2)- 
Finally, the ratio of the vortex core radius to the intervortex separation follows from 

and use of Eq. (|2^) to eliminate 77 in favor of Q gives the generalization of Eq. ([l4]) for the case of a combined quadratic 
and quartic confining potential. Figure || shows the f2 dependence of £, 2 /b 2 for two illustrative cases; the sharp increase 
near f2 C 2 and the finite limiting value are particularly evident. 



V. DISCUSSION AND CONCLUSIONS 



The behavior of a dense vortex lattice in a dilute trapped rotating Bose-Einstein condensate has acquired a new 
interest because of recent observations of such lattices The experiments use radial harmonic traps, and the loss 

of confinement as the external rotation speed f2 approaches the trap frequency to± implies that the behavior becomes 
singular As a result, it would be difficult to study the approach to what is effectively the upper critical angular 
velocity Vt C 2 ~ fr/Mt; 2 when the vortex cores overlap. The present analysis considers the addition of a stiffer quartic 
potential V±(r±) = \kr^, which ensures confinement for any O. Thus, the approach to f2 C 2 is more gradual and 
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could well be observed. In particular, f2 C 2/w_L exceeds 1 and should be accessible to experiments. In principle, a 
similar analysis is possible for other stiff radial confining potentials, for example, other power laws or an optical dipole 
waveguide made from a hollow blue-detuned laser beam ]^3| . 

The Thomas-Fermi approximation assumes that the total kinetic energy of the condensate is much smaller than the 
energies associated with the external trapping potential and the interparticle Hartree potential. This picture certainly 
applies for a relatively small number of vortices when £ <C b <C R±, but it fails near J7 C 2 when £ < b because of the 
density gradient near the many vortex cores. This effect will not qualitatively alter the value of f2 C 2 determined for 
the quartic potential, but it will affect the detailed description for f2 < Q C 2- This interesting question remains open. 

Another question concerns the number of atoms No in the condensate, which here has been assumed to remain 
fixed. This picture also must fail sufficiently near f2 C 2 because the vortex cores fill the entire volume. For type-II 
superconductors Q| at a given temperature, the vortex core size £ remains fixed. In that system, an increased applied 
magnetic field reduces the intervortex separation b, leading to the disappearance of the superconducting component 
because the "normal" cores eventually overlap. In a dilute Bose-Einstein gas, however, the core size £ = h/y/2Mfi 
itself increases and diverges as fi — > tt C 2 because of the decreasing chemical potential. Thus the approach to the upper 
critical angular velocity is more sudden in a dilute trapped gas. Unfortunately, inclusion of the quantum depletion is 
fairly complicated for a spatially nonuniform medium pij-pq] although the general formalism is well known p9| , pp| . 

At zero temperature in the grand canonical ensemble for a rotating system in equilibrium at chemical potential fi 
and an angular velocity f2, for example, the ground-state expectation value of the operator K = H — fiN — QL Z is the 
thermodynamic potential F(fj,,il) — (K). In the presence of a Bose-Einstein condensate with Ao condensate atoms, 
this function also depends on Ao, so that (K) = F(/x, fi; A^)- The usual thermodynamic relation A(/z, Q; No) = 
— (dF/d[i)n]\[ determines the mean number of particles, which here depends not only on /j, and fi, but also on Ao. 
This latter parameter can be fixed by adjusting Ao to minimize F, so that {dF / dNo) = 0. For a uniform Bose gas 
in a stationary box of volume V , it is straightforward to verify that this procedure yields the familiar zero-temperature 
depletion of the condensate (N — No)/N w ^{na 3 /it) 1 / 2 as well as the first correction to the chemical potential [^9] [[l). 
It should be feasible to extend this analysis to a dilute rotating trapped Bose condensate, and this problem definitely 
merits further study. 
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FIG. 1. Dimensionless chemical potential Jl = n/hu}± for the combined quadratic and quartic radial confining potential as 
a function of the dimensionless angular velocity f2 = ft/ui± with (Noa/d±) (lo z /lj±) — 10 4 . The two curves correspond to the 
values (a) di/d± — 5 (b) di/d± — 2. Note that Jl remains finite at Q = 1 for both values of d^jd^. 



FIG. 2. Dimensionless squared radius R\jd\ for the combined quadratic and quartic radial confining potential as a function 
of the dimensionless angular velocity Q = Q/lo± with (Noa/d±) (lo z /uj±) = 10 4 . The two curves correspond to the values (a) 
d4,/d± = 5 (b) d4/d± — 2. In contrast to the behavior seen in Eq. ( |lo| ) for a rotating harmonic radial potential, the squared 
radius here remains finite for both values of d4,/d±_ as Q, — > S7 C 2 given in Eq. (0). 



FIG. 3. Fraction of volume (, 2 /b 2 occupied by the vortex cores for the combined quadratic and quartic radial confining 
potential as a function of the dimensionless angular velocity Q = Q,/ljj_ with (Noa/d±_) (lu z /lo±) = 10 4 . The two curves 
correspond to the values (a) cLi/dx = 5 (b) d4/d± — 2. Note that £ 2 /b 2 remains finite for both values of di/dx as Q — > Q C 2 
given in Eq. (p^). 
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